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1.  INTRODUCTION 


While  finite  element  hydrocodes  can  predict  the  contours  of  finitely  deforming  metals  with 
reasonable  accuracy,  their  ability  to  determine  local  strain  and  stress  states  leading  to 
catastrophic  failure  by  such  mechanisms  as  adiabatic  shear  banding  and  void  openings  leaves 
much  to  be  desired.  An  improvement  in  the  ability  to  compute  such  local  state  histories  would 
significantly  enhance  the  design  of  warheads  and  predictions  of  armor  penetration  and 
behind-armor  debris  formation.  In  recognition  that  better  material  characterization  was  an 
essential  ingredient  of  the  desired  improvement,  a  joint  U.S.  Army  Ballistic  Research 
Laboratory/Materials  Technology  Laboratory  (BRL/MTL)  program  entitled  Advanced 
Constitutive  Models  was  initiated  several  years  ago.  This  report  discusses  one  facet  of  this 
investigation. 

At  BRL,  it  was  decided  to  employ  the  DYNA3D  hydrocode  (Hallquist  1983)  in  connection 
with  this  study.  This  is  a  Lagrangian  finite  element  computer  program  currently  in  widespread 
use  which  employs  an  adequate  finite  deformation  formulation  and  features  a  choice  of 
approximately  30  constitutive  models.*  However,  only  a  few  of  these  models  are  suitable  for 
the  large  strain,  low  rate  applications  to  be  discussed  in  sequel.  These  models  generally  ** 
employ  the  von  Mises  yield  condition  Jg  -  =  0,  where  Jg  is  the  second  invariant  of  the 

stress  deviator  tensor,  as  the  plastic  potential  function  and  use  the  Jaumann  stress  rate  to 
account  for  the  rigid  body  rotation  of  elements. 

2.  QUASI-STATIC  MATERIAL  PROPERTIES 

For  an  evaluation  of  material  constitutive  behavior,  it  is  necessary  to  have  test  data  for  the 
specific  lot  of  material  to  be  characterized.  It  was  decided  that  this  investigation  would 
commence  with  a  study  of  oxygen-free,  high  conductivity  (OFHC)  copper  and  MTL  was  tasked 


*  This  term  is  used  to  identify  a  mathematical  function  or  algorithm  which  determines  stress  tensor  components  from 
the  current  rate-of-deformation  components,  some  past  history  data,  and  possibly  temperature  through  an  incremental 
plasticity  time  marching  process. 

**  An  exception  is  the  model  of  D.  Bammann  (Johnson  and  Bammann  1984)  as  incorporated  by  M.  Chiesa  of  Sandia/ 
Livermore  which  is  a  multi-parameter  micromechanically-based  constitutive  system  which  uses  the  Green-Naghdi  stress 
rate. 
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to  perform  the  necessary  experimentation.  Dr.  Tusit  Weerasooriya  has  reported 
(Weerasooriya  and  Swanson  1991)  data  from  two  types  of  quasi-static  (isothermal)  tests  on 
annealed  copper: 

(1)  Uniaxial  compression  stress-strain  curves  for  the  range  -1.30  <  <  0  natural  strain. 

(2)  Torsional  shear  stress-strain  data  from  twist  tests  on  modified  Lindholm-type 

thin-walled  tubular  specimens  for  the  range  0  <  <  1.4  (tensor)  shear  strain.  These 

tests  were  performed  for  two  conditions  of  axial  restraint:  (a)  almost  total  axial 
restraint  in  which  case  the  induced  axial  force  was  recorded  and  (b)  no  applied  axial 
restraint  where  the  axial  displacement  was  monitored. 

Also,  Dr.  Weerasooriya  provided  the  author  a  curve  for  reversed  loading  of  a  torsion  specimen 
which  permitted  an  assessment  of  the  Bauschinger  effect  for  this  material. 

3.  CONSTITUTIVE  MODELING  CONSIDERATIONS 

The  compression  test  data  cited  above  reveal  that  the  stress-strain  curve  for  annealed 
copper  is  nonlinear  over  the  entire  range,  the  elastic  portion  being  of  negligible  size.  For  a 
curve  of  this  form,  a  bilinear  representation  (such  as  DYNA3D  Material  3)  is  unsatisfactory. 
Material  10  of  the  DYNA  code,  which  permits  input  of  up  to  16  stress-strain  points  and 
interpolates  linearly  for  intermediate  values,  is  more  appropriate  but  only  treats  isotropic 
hardening. 

The  Lindholm-type  torsion  specimen  does  not  result  in  a  homogeneous  state  of  stress  in 
the  thin-walled  test  section;  also,  the  shearing  strain  and  plastic  deformation  extend  into  the 
transition  section.  For  this  reason,  various  investigators  have  resorted  to  3-D  finite  element 
mode  ng  of  the  whole  specimen  (using  an  assumed  constitutive  model)  to  provide  a  basis  for 
interpretation  of  test  data.  Lipkin  et  al.  (1987)  reported  a  DYNA3D  calculation  using  an  earlier 
version  of  the  Bammann  constitutive  model  in  which  a  twist  rate  high  enough  to  predict 
adiabatic  shear  banding  was  used.  The  present  author  performed  a  DYNA3D  analysis  for  the 
MTL  geometry  using  Material  3  and  found  this  geometry  was  prone  to  premature  torsional 
buckling.  It  was  recommended  that  a  thicker  wall  be  used  since  an  interpretive  analysis  would 
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be  required  in  any  event.  Dr.  C.  S.  White  (1990)  reported  an  ABAQUS  anaiysis  of  a 
geometric  configuration  cioseiy  corresponding  to  that  used  by  Weerasooriya  and  has 
concluded  that  "about  78%  of  the  twist  that  is  applied  at  the  grips  actually  goes  into  the 
deformation  in  the  gauge  section."  Although  White’s  calculations  were  made  for  a  different 
material,  it  was  decided  to  multiply  Weerasooriya’s  shear  strain  data  by  a  factor  of  0.78. 

4.  USE  OF  EXPERIMENTAL  DATA;  DISCREPANCIES 

If  the  compression  test  data  are  used  as  input  to  DYNA3D  Material  1 0  and  a  calculation  is 
made  for  an  element  constrained  to  deform  uniaxially  with  no  lateral  restraints,  the  code 
predictions  are  (naturally)  in  excellent  agreement  with  the  experimental  values.  (It  is 
necessary  to  choose  an  equation-of-state  which  permits  specifying  the  pressure  to  be 
proportional  to  the  volumetric  strain  (proportionality  constant  =  bulk  modulus),  use  the 
hourglass  viscosity  type  3  {Flanagan-Belytschko  (1981)  with  exact  volume  integration),  and  to 
set  the  Gruneisen  coefficient  =  0  for  an  isothermal  calculation.)  Similarly,  when  torsion  test 
data  are  inserted  in  DYNA3D  Material  10  and  a  simple  shear  problem  is  run  for  an  element, 
there  is  no  discrepancy. 

However,  when  compression  test  data  are  used  in  Material  1 0  to  predict  the  stresses  in  an 
element  subjected  to  increasing  simple  shear  the  result  shown  in  Figure  1  is  obtained,  where 
the  overprediction  of  the  shear  stress  is  as  great  as  30%.  Not  surprisingly,  the  converse  is 
also  true:  use  of  the  torsion  test  data  as  input  for  a  uniaxial  compression  calculation  results  in 
a  significant  underprediction  of  the  axial  stress  as  shown  in  Figure  2. 

This  phenomenon  is  no  new  discovery  and  has  been  discussed  in  the  literature  by  Prager 
(1945),  Drucker  (1949),  Edelman  and  Drucker  (1951),  and  many  others.  The  basic  problem  is 
that  use  of  the  von  Mises  condition  as  a  loading  function  for  work  hardening  materials  does 
not  closely  describe  the  plastic  deformation  of  many  materials  (even  though  it  is  widely 
employed  for  this  purpose  in  most  currently  used  hydrocodes).  According  to  the  authors  just 
cited,  the  problem  can  be  resolved  by  use  of  a  loading  function  depending  on  both  the  second 
and  third  invariants  of  the  stress  deviation;  i.e.,  a  (Jg,  Jg)  theory.  Although  such  theories  seem 
to  have  fallen  into  disuse,  perhaps  this  approach  should  be  considered  for  applications 
involving  large  strains,  where  the  discrepancies  are  greatest. 
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Ea  (TMUI  SnUM) 

P'Qure  1 .  Simple  Shear  Calculation  (Geometric  Constraints.  Material  Tvoe  1QV 


Uniaxial  Compression  Calculation  (Geometric  Constraints.  Material  Tvoe  1Q^ 
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5.  IMPLEMENTATION  OF  A  (Jg.  J3)  THEORY 


It  was  decided  to  employ  the  function  quoted  by  Malvern  (1969): 


f(J2,J3)*J; 


1  - 


C(J3)^ 

(Ja)® 


-  =  0 


where 


(1) 


(2) 


'Ja  *  S,,  S22  833  +  2S,2  S23  S3,  -  S,,  823  -  S22  S31  -  S33  S12 


=  8,  S2  S3 


(3) 


k  =  Yield  stress  in  pure  shear,  variable  for  a  work  hardening  material 

c  =  nondimensional  parameter  to  be  adjusted  to  provide  match  of  post  yield 
flow  data 

Sjj  =  deviatoric  stress  components 

% 

8^  =  principal  deviatoric  stresses 


Installation  of  Equation  1  as  the  load  function  in  DYNA3D  was  easily  accomplished  since 
current  values  of  ail  required  quantities  are  available  in  the  stress  evaluation  subroutine. 
Because  J3  vanishes  for  pure  shear,  it  is  attractive  to  insert  the  tabular  data  from  the  shear 
tests  in  Material  1 0  and  then  determine  the  parameter  c  to  provide  correspondence  to  the 
compression  test  data  for  uniaxial  stress  calculations.  For  the  present  application,  the  peak 
stress  was  employed  because  the  uniaxial  stress-strain  curve  appears  to  level  off  at  this  value 
and  the  interest  is  in  large  strain  plasticity.  Of  course,  other  matching  criteria  could  be 
adopted  for  other  strain  ranges.  Only  a  few  trials  were  required  to  obtain  c  =  2.64;  this  value 
was  used  in  DYNA3D  to  obtain  the  solid  curves  displayed  in  Figures  1  and  2.  It  may  be  seen 
in  Figure  2  that  there  is  a  significant  difference  between  experimental  and  predicted  values  in 
the  small  strain  region.  This  is  attributed  to  the  much  more  rounded  "knee"  exhibited  by 
copper  in  compression  than  in  tension  or  shear. 
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6.  GEOMETRIC  REPRESENTATION 


It  is  instructive  to  view  the  (Jg,  J3)  yield  function  in  principal  stress  space,  where  it  may  be 
recalled  that  the  von  Mises  function  is  represented  by  a  circular  cylinder  coaxial  with  the 
hydrostatic  line  0^  =  02  =  03.  To  achieve  this,  the  following  orthogonal  coordinate 
transformation  was  made; 


Oi  = 


s/S  s/2 


O2 


O3 


\l2  sjz 

^  T3 

/6  s[Z 


(4) 


With  this  transformation,  the  T3  axis  coincides  with  the  hydrostatic  line  and  the  t,  axis  is 
aligned  with  the  projection  of  the  03  axis  on  the  deviatoric  plane  a,  +  Og  +  03  =  0.  Using 
Equations  4,  one  obtains 


^^2  (^1 

(5) 

% 

■  t,T|) 

(6) 

v/6  3 

and,  substituting  these  values  in  Equation  1 ,  there  results 


27X2  +  {(81  -36c)x?-54k2)x2  +  {(81  +  2Ac)xl 
+  (27-4c)x?-54k2xt  =  0 


(7) 


which  defines  the  contour  of  the  yield  function  in  the  x,.  Xg-plane. 

As  noted  by  Hill  (1950),  it  is  only  necessary  to  compute  coordinates  for  a  30°  segment  of 
this  plane,  the  remainder  of  the  locus  being  determined  by  the  symmetry  constraints  for  an 
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isotropic  material.  This  locus  for  c  =  2.64  is  displayed  in  Figure  3  as  a  solid  line.  It  is  seen 
that  the  (J2,  J3)  surface  is  a  fluted  cylinder  with  the  von  Mises  cylinder  inscribed.  For  Material 
1 0  which  provides  only  for  isotropic  hardening,  these  surfaces  would  expand  uniformly  as 
plastic  deformation  progresses.  It  should  be  mentioned  that  DYNA3D  presently  uses  only  the 
Krieg-Key  radial  return  algorithm  (Krieg  and  Key,  1976)  to  return  the  stress  state  to  the 
updated  yield  surface  each  cycle.  This  algorithm  is  clearly  more  appropriate  for  the  Jg  surface 
than  for  the  (J2,  J3)  function.  However,  it  is  still  attractive  due  to  its  simplicity  and  the 
assurance  that  the  yield  surface  will  be  intersected,  although  for  the  latter  it  corresponds  to  a 
nonassociated  flow  rule.*  A  "normal  return"  algorithm  would  seem  preferable  or  perhaps  the 
recently  published  Nemat-Nasser  algorithm  (Nemat-Nasser  1991),  but  either  of  these  would 
require  more  extensive  calculations  per  cycle.  These  matters  may  be  academic  in  view  of  the 
following  discussion. 

7.  MATERIAL  STABILITY,  CONVEXITY  REQUIREMENTS 

It  is  clear  from  inspection  that  the  c  ®  2.64  loading  function  in  Figure  3  violates  the 
requirement  that  the  surface  be  convex,  as  deduced  by  Drucker  (1951,  1959)  from  his 
postulates  for  material  stability.  Although  the  author  has  made  a  considerable  number  of 
computer  runs  using  this  value  of  c  with  no  evidence  of  instability,  these  cases  were  not 
designed  to  test  all  types  of  loading  and  unloading.  Thus,  it  is  appropriate  to  ask:  Can  any 
(J2,  J3)  load  function  satisfy  the  convexity  requirement?  The  answer  is  yes  and  the  limitation 
on  c  is  determined  by  evaluating  the  curvature  of  the  load  function  at  its  point  of  tangency  to 
the  von  Mises  circle.  Thus,  at  =  0, 


dx?  ^2  k 


(8) 


This  happens  to  be  the  exact  curvature  since _ i  =  0  at  this  location. 

dx, 


*  This  consideration  does  not  invalidate  the  results  displayed  in  this  report  since,  for  both  pure  shear  and  uniaxial 
loading,  radial  and  normal  return  coincide. 
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) 


Figure  3.  Axial  View  of  Yield  Loci. 

Therefore,  the  transition  in  sign  of  the  curvature  occurs  at  c  »  0.75  and  the  load  function  is 
everywhere  convex  for  this  or  lesser  values  of  c.  The  form  of  the  load  function  for  c  =  0.75  is 
plotted  in  the  upper  sextant  of  Figure  3.  Unfortunately,  for  this  value  of  c,  the  percentage 
difference  between  the  maximum  and  minimum  radii  of  the  load  function  is  only  about  6%.  It 
can  now  be  seen  that  the  greatest  percentage  difference  in  radii  consistent  with  convexity 
requires  a  load  function  in  the  form  of  a  hexagon  circumscribed  about  the  von  Mises  circle,  as 
shown  in  the  upper  portion  of  Figure  3.  Even  for  this  form  of  yield  function  (with  its  attendant 
analytical  complexities),  the  percentage  difference  in  radii  is  only  1 5.5%,  about  half  of  what  is 
needed  to  reconcile  the  discrepancies  associated  with  experimental  data. 

8.  CONCLUDING  REMARKS 

It  is  realized  that  other  investigators  (Shrivastava,  Jonas,  and  Canova  1982;  Asaro  and 
Needleman  1985;  Weerasooriya  and  Swanson  1991)  have  treated  the  subject  matter  of  this 
paper  by  micromechanical  modeling  and  have  attributed  the  cited  discrepancies  to  formation 
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of  texture.  However,  the  end  objective  of  such  research  does  not  appear  to  be  the 
identification  of  load  functions  appropriate  to  classical  plasticity. 

This  author  is  convinced  that  material  behavior  must  be  modeled  within  the  framework  of 
continuum  mechanics  if  efficient,  large-scale  computations  are  to  yield  valid  results  for 
engineering  purposes.  While  the  study  reported  herein  may  seem  inconclusive,  it  is  believed 
that  use  of  a  nonconvex  loading  function  may  be  permissible  for  certain  classes  of  problems, 
especially  if  the  computer  program  is  modified  to  test  on  the  sign  of  the  plastic  work  increment 
in  each  element  and  to  terminate  the  calculations  if  a  negative  work  increment  is  predicted.  It 
is  very  desirable  to  resolve  this  matter  before  proceeding  to  studies  of  strain-rate  effects, 
stress  rate  models,  etc.,  where  uncertainties  regarding  the  loading  function  may  obscure 
interpretation  of  other  types  of  experimental  results. 

It  should  be  mentioned  that  the  author  has  added  a  new  subroutine  to  a  research  version 
of  DYNA3D  which  includes  many  effects  previously  not  present  in  a  single  material  model, 
specifically: 

(1)  Finite  elastoplastic  straining. 

% 

(2)  Mixed  kinematic/isotropic  hardening,  as  recommended  by  Hodge  (1957). 

(3)  Arbitrary  shape  of  uniaxial  stress-strain  function  through  input  of  tabular  data. 

(4)  Choice  of  stress  rate  formulation  (Jaumann  [1 905]  or  Green  and  Naghdi  [Green  and 
Naghdi  1965;  Green  and  Mclnnis  1967]). 

(5)  Choice  of  several  equation-of-state  models. 

It  is  planned  to  use  this  tool  to  study  finite  plasticity  for  various  proportional  and 
nonproportional  loading  paths. 


9 


Intentionally  left  blank. 


10 


9.  REFERENCES 


Asaro,  R.  J.,  and  A.  Needleman.  "Texture  Development  and  Strain  Hardening  in  Rate 
Dependent  Polycrystals."  Acta  Metall..  vol.  33,  p.  923,  1985. 

Drucker,  D.  C.  "Relation  of  Experiments  to  Mathematical  Theories  of  Plasticity."  Journal  of 
Applied  Mechanics,  vol.  16,  no.  4,  pp.  349-357,  1949. 

Drucker,  D.  C.  "A  More  Fundamental  Approach  to  Plastic  Stress-Strain  Relations."  Proc.  First 
U.S.  Natl.  Cong.  Appl.  Mechanics,  ASME,  pp.  487-491,  1951. 

Drucker,  D.  C.  "A  Definition  of  Stable  Inelastic  Material."  Journal  of  Applied  Mechanics. 
vol.  26.  pp.  101-106,  1959. 

Edelman,  F.,  and  D.  C.  Drucker.  "Some  Extensions  of  Elementary  Plasticity  Theory."  Journal 
of  the  Franklin  Institute,  vol.  251,  pp.  581-605,  1951. 

Flanagan,  D.  P.,  and  T.  Belytschko.  "A  Uniform  Strain  Hexahedron  and  Quadrilateral  and 
Orthogonal  Hourglass  Control."  Int.  J.  Numer.  Meths.  Eno..  vol.  17,  pp.  679-706,  1981. 

Green,  A.  E.  and  P.  M.  Naghdi.  "A  General  Theory  of  an  Elastic-Plastic  Continuum."  Arch. 
Rat.  Mech.  Anal.,  vol.  18.  pp.  251-281,  1965. 

Green,  A.  E.  and  B.  C.  Mclnnis.  "Generalized  Hypo-Elasticity."  Proc.  Roy.  Soc.  Edinburgh. 
A57.  p.  220,  1967. 

Hallquist,  J.  O.  "Theoretical  Manual  for  DYNA3D."  UCID-19401,  University  of  California, 
Lawrence  Livermore  National  Laboratory,  1983. 

Hill.  R.,  "The  Mathematical  Theory  of  Plasticity."  Oxford,  p.  18,  1950. 

Hodge,  P.  G.,  Jr.  Discussion  of  Prager  (1956),  Journal  of  Applied  Mechanics,  vol.  24,  no.  3, 
pp.  482-483,  1957. 

Jaumann,  G.  "Grundlagen  der  Bewegungslehre."  Leipzig,  1905. 

Johnson,  G.  C.,  and  D.  J.  Bammann.  "A  Discussion  of  Stress  Rates  in  Finite  Deformation 
Problems."  Int.  J.  Solids  Structures,  vol.  20,  no.  8,  pp.  725-737,  1 984. 

Krieg,  R.  D.  and  S.  W.  Key.  "Implementation  of  a  Time  Independent  Plasticity  Theory  Into 
Structural  Computer  Programs."  Constitutive  Equations  in  Viscoplasticity,  ASME,  AMD 
vol.  20,  1976. 

Lipkin,  J.,  M.  L.  Chiesa,  and  D.  J.  Bammann.  "Thermal  Softening  of  304L  Stainless  Steel: 
Experimental  Results  and  Numerical  Simulations."  Proc.  of  IMPACT’87,  Bremen, 

Germany,  1987. 


11 


Malvern,  L.  E.  “Introduction  to  the  Mechanics  of  a  Continuous  Medium."  Prentice-Hall, 
p.  355,  1969. 

Nemat-Nasser,  S.  "Rate-Independent  Finite-Deformation  Elastoplasticity:  A  New  Explicit 
Constitutive  Algorithm.  Mechanics  of  Materials,  vol.  11,  No.  3,  pp.  235-249,  1991. 

Prager,  W.  "Strain  Hardening  Under  Combined  Stresses."  J.  AddI.  Physics,  vol.  1 6, 
pp.  837-640,  1945. 

Shrivastava,  S.  C.,  J.  J.  Jonas,  and  G.  Canova.  "Equivalent  Strain  in  Large  Deformation 
Torsion  Testing:  Theoretical  and  Practical  Considerations."  J.  Mech.  Phvs.  Solids. 
vol.  30.  pp.  75-90,  1982. 

Weerasooriya,  T.,  and  R.  A.  Swanson.  "Experimental  Evaluation  of  the  Taylor-Type 
Polycrystal  Model  for  the  Finite  Deformation  of  an  FCC  Metal  (OFHC  Copper)." 

MTL  TR  91-20,  U.S.  Army  Materials  Technology  Laboratory,  1991. 

White.  C.  S.  "An  Analysis  of  the  Torsion  Specimen  Used  in  Constitutive  Modeling." 
Presented  at  the  Eighth  Army  Conference  on  Applied  Mathematics  and  Computing, 
Cornell  University,  1990. 


12 


No.  of 

Copies  Organization 


No.  of 

Copies  Organization 


2  Administrator 

Defense  Technicai  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22304-6145 


1  Commander 

U.S.  Army  Tank- Automotive  Command 
ATTN;  ASQNC-TAC-DIT  (Technical 
Information  Center) 

Warren,  Ml  48397-5000 


1  Commander  1 

U.S.  Army  Materiel  Command 
ATTN:  AMCAM 
5001  Eisenhower  Ave. 

Alexandria,  VA  22333-0001 

1 

1  Commander 

U.S.  Army  Laboratory  Command 
ATTN:  AMSLC-DL 
2800  Powder  Mill  Rd. 

AdelphI,  MD  20783-1145  ((Waniyn 

2  Commander 

U.S.  Army  Armament  Research, 

Development,  and  Engineering  Center 
ATTN:  SMCAR-IMI-I  (Und-..  oniy)i 

Picatinny  Arsenal,  NJ  07806-5000 

2  Commander 

U.S.  Army  Armament  Research, 

Development,  and  Engineering  Center  1 

ATTN:  SMCAR-TDC 
Picatinny  Arsenal,  NJ  07806-5000 


1  Director 

Benet  Weapons  Laboratory  2 

U.S.  Army  Armament  Research, 

Development,  and  Engineering  Center 
ATTN:  SMCAR-CCB-TL 

Watervliet,  NY  12189-4050  l 

(UnctaM.  ooiy)i  Commander 

U.S.  Army  Rock  Island  Arsenal  3 

ATTN;  ^CRI-TL/Technical  Library 
Rock  Island,  IL  61299-5000 

1  Director 

U.S.  Army  Aviation  Research  1 

and  Technology  Activity 
ATTN:  SAVRT-R  (Library) 

M/S  219-3  10 

Ames  Research  Center 
Moffett  Field,  CA  94035-1000 


Director 

U.S.  Army  TRADOC  Analysis  Command 
ATTN"  ATRC-WSR 

White  Sands  Missile  Range,  NM  88002-5502 
Commandant 

U.S.  Army  Field  Artillery  School 

ATTN:  ATSF-CSI 

Ft.  Sill,  OK  73503-5000 

Commandant 

U.S.  Army  Infantry  School 

ATTN:  ATSH-CD  (Security  Mgr.) 

Fort  Benning,  GA  31905-5660 

Commandant 
U.S.  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  QA  31905-5660 

WL/MNOI 

Eglin  AFB,  FL  32542-5000 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 

Cdr,  USATECOM 
ATTN:  AMSTE-TC 

Cdr,  CRDEC,  AMCCOM 
ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-MSI 

Dir,  VLAMO 
ATTN:  AMSLC-VL-D 

Dir,  USABRL 
ATTN:  SLCBR-DD-T 


1  Commander 

U.S.  Army  Missile  Command 
ATTN:  AMSMI-RD-CS-R  (DOC) 
Redstone  Arsenal,  AL  35898-5010 


13 


No.  of 


No.  of 

Coptes  Organization 

2  Director 
DARPA 

ATTN:  J.  Richardson 
LTC  J.  Beno 
1400  Wilson  Blvd. 

Arlington,  VA  22209-2308 

1  Defense  Nuclear  Agency 
ATTN:  MAJ  James  Lyon 
6801  Telegraph  Road 
Alexandria,  VA  22192 

4  U.S.  Army  Research  Office 
ATTN:  I.  Ahmad 

K.  Iyer 
J.  Wu 

Technical  Library 
P.O.  Box  12211 
4300  Miami  Blvd. 

Research  Triangle  Park,  NC  27709 

3  U.S.  Army  MICOM 
ATTN:  AMSMI-RD-ST-WF, 

Lynn  Craft 
Donald  Lovelace 
Michael  Schexnayder 
Redstone  Arsenal,  AL  35898-5250 

2  U.S.  Army  Belvoir  RD&E  Center 
ATTN:  S.  G.  Bishop 

C.  Kominos 

Fort  Belvoir,  VA  22060-5166 

5  Commander 

U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  T.  Davidson 

V.  Lindner 
J.  Pearson 
E.  Baker 
Technical  Library 

Picatinny  Arsenal,  NJ  07806-5000 


Copies  Oroangation 

6  U.S.  Army  Materials  Technology  Laboratory 
ATTN:  SLCMT-MRT, 

S. C.  Chou 

C. S.  White 

J.  McLaughlin 

T.  Weerasooriya 

D.  Dandekar 

Technical  Lbrary 

Watertown,  MA  02172-0001 

2  Naval  Weapons  Center 

ATTN:  Don  Thompson,  Code  3268 
Technical  Library 
China  Lake,  CA  93555 

2  Commander 

Naval  Surface  Warfare  Center 
ATTN:  William  Mock 

Technical  Library 
Dahlgren,  VA  22448-5000 

7  Commander 

Naval  Surface  Warfare  Center 
ATTN:  P.  C.  Huang,  G-402 
J.  P.  Matra 
P.  Walter 
L.  Mensi 
F.  J.  Zerilli 
H.  Chen,  U-12 
Technical  Library 
10901  New  Hampshire  Ave. 

Silver  Spring,  MD  20903-5000 

1  Sandia  National  Laboratories 
ATTN:  D.  Bammann 
J.  Lipkin 

Livermore,  CA  94550 

1  Air  Force  Wright  Laboratories 
Materials  Laboratory 
ATTN:  Dr.  T.  Nicholas 
Wright-Patterson  AFB,  OH  45433 

3  Air  Force  Armament  Laboratory 
ATTN:  AFATUDUW, 

W.  Cook 
AFATL/MNW, 

J.  Foster 
Technical  Library 
Eglin  AFB,  FL  32542 


14 


Director 

SarKlia  National  Latx)ratories 
ATTN:  J.M.  McGlaun 
P.  Yarrington 
E.  Hertel 
Technical  Library 
P.O.  Box  5800 
Albuquerque,  NM  87185 

Director 

Los  Alamos  National  Laboratory 
ATTN:  G.E.  Cort,  F663 

T. F.  Adams,  F663 
D.  Mandell,  F663 
R.  Karpp,  MS  J960 
J.  Chapyak 
Technical  Library 

P.O.  Box  1663 

Los  Alamos,  NM  87545 

Director 

Lawrence  Livermore  National  Laboratory 
ATTN:  R.  Tipton.  L-35 
R.  Whirley.  L-122 
R.  Christensen,  L-35 
Technical  Library 
P.O.  Box  808 
Livermore.  CA  94550 

Aerojet  Precision  Weapons 
Dept.  5131/T-W 
ATTN:  J.  Carleone 
1100  Hollyvale 
Azusa,  CA  91702 

Dyna  East  Corporation 
ATTN:  P.C.  Chou 
R.  Ciccarelli 
W.  Flis 

3201  Arch  Street 
Philadelphia,  PA  19104 

Southwest  Research  Institute 
ATTN:  C.  Anderson 
A.  Wenzel 

U.  Lindholm 
P.O.  Drawer  28255 

San  Antonio,  TX  78228-0255 


1  Alliant  Techsystems,  Inc. 

ATTN:  G.R.  Johnson 
MN  48-2700 
7225  Northland  Drive 
Brooklyn  Park,  MN  55428 

1  S-Cubed 

ATTN:  R.  Sedgwick 

P.O.  Box  1620 

La  Jolla.  CA  92038-1620 

2  Orlando  Technology,  Inc. 

ATTN:  D.  Matuska 

J.  Osborn 
P.O.  Box  855 
Shalimar,  FL  32579 

1  Livermore  Software  Technology  Corp. 

ATTN:  John  O.  Hallquist 
2876  Waveriy  Way 
Livermore,  CA  94550 

1  Rensselaer  Polytechnic  Institute 

ATTN:  Professor  E.  KrempI 
Troy,  NY  12181 

3  Brown  University 
Division  of  Engineering 
ATTN:  Professor  R.  Clifton 

Professor  H.  Kolsky 
Professor  R.  Asaro 
Providence,  Rl  02912 

1  Camegie-Mellon  University 
Department  of  Mathematics 
ATTN:  Dr.  M.E.  Gurtin 
Pittsburgh,  PA  15213 

2  The  Johns  Hopkins  University 
ATTN:  Professor  W.  Sharpe 

Professor  J.  F.  Bell 
34th  and  Charles  Streets 
Baltimore,  MD  21218 

1  University  of  California  at  San  Diego 

Department  of  Mechanical  and  Aerospace 
Engineering 

ATTN:  Professor  S.  Nemat-Nasser 
La  Jolla,  CA  92093 


No.  of 

Copies  Oroanization 

1  University  of  Delaware 

Department  of  Mechanical  and  Aerospace 
Engineering 

ATTN:  Professor  J.  Vinson 
Newark.  DE  19711 

2  University  of  Florida 

Department  of  Engineering  Science  and 
Mechanics 

ATTN;  Professor  L.  Malvern 
Professor  D.  Drucker 
Gainesville,  FL  32601 

1  Virginia  Polytechnic  Institute  and  State 
University 

Department  of  Engineering  Science  and 
Mechanics 

ATTN:  Professor  C.  W.  Smith,  Jr. 
Blacksburg,  VA  24061 


16 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it 
publishes.  Your  comments/answers  below  will  aid  us  in  our  efforts. 

1.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of 
interest  for  which  the  report  will  be  used.)  _ 


2.  How,  ^ecifically,  is  the  report  being  used?  (Information  source,  design  data,  procedure, 
source  of  ideas,  etc.)  _ 


3.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or 
dollars  saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please 
elaborate.  _ 


4.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports? 
(Indicate  changes  to  organization,  technical  content,  format,  etc.)  _ 


BRL  Report  Number  BRL-tr-3368 _  Division  Symbol 

Check  here  if  desire  to  be  removed  from  distribution  list.  _ 

Check  here  for  address  change.  _ 

Current  address:  Organization  _ 

Address  _ 


Department  of  the  army 

Director 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-5066 


OFRCIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

RRST  CLASS  FERMI  No  0001,  APG.  MO 


Postage  wItS  be  paid  by  addressee 


NO  POSTAGE 
NECESSARY 
IF  MAILEO 
IN  THE 

UNITED  STATES 


Director 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground.  MD  21005-5066 


